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The  steady  two-dlnen«ional  sink '•type  flow  of  a viscous^  heat- 
conducting.  perfect  gHa  i3  invc  atigated , Axi  approximate  solution  of  this 
prrbl.'-'n  i*  given  for  the  case  ol  Itrge  Reynolds  number  He  (c£.  ti»e 
definition  g:iven  in  the  text).  In  obtaining  the  present  &t'.'.uilinn  the  values 
of  Prardtl  ntimber  and  the  ratio  of  the  first  and  second  viecocity 
coefficient  may  be  arbitrary.  The  result  shows  that  the  solution  has  two 
branches,  both  of  physical  significance.  On  the  subsonic  branch  of  the 
solution  the  flew  speed  starts  from  ilie  stagnation  point  at  infinity,  increases 
monotonically  for  decreasing  radial  distance  and  everraally  terminates 
with  maximum  speed  at  a certain  point  inside  the  inviscid  sonic  circle. 

The  solutions  of  the  supersonic  branch,  which  start  with  the  maximum 

speed  at  infinity,  all  coiitain  cylindrical  shocks.  Within  the  shock  the 

flow  speed  assumes  a minimum  value  and  after  the  shock  all  solutions 

tend  asymptotically  to  the  subsonic  br&nch.  In  contrast  to  the  plane  shock, 

the  cylindrical  shock  strength  is  limited  to  the  order  0( 

the  shock -thickness,  of  . The  latter  quantity  implies  that  the 

thickness  of  the  region  in  which  the  viscous  effects  are  important  is 

thinner,  in  order  of  magnitude,  than  that  of  ordinary  boxindary  layer 

(of  )•  but  is  thicker  than  that  of  plane  normal  shock  (of 

It  is  found  that  the  entropy  of  the  supersonic  branch  rises  to  a maximum 

within  the  shock  while  for  the  subsonic  branch,  the  entropy  increases 


monotonically  with  the  radial  distance.  The  total  variation  of  the  entropy 
across  the  shock  is  found  to  be  of  > which  is  in  general  greater 

than  that  across  a plane  normal  shock  Q (shock  strength^  J . The 
effect  on  the  flow  quantities  due  to  the  variation  in  viscosity  coefficients, 
assumed  to  depend  on  the  local  temperature,  is  found  to  be  at  most  of 
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Introduction 


I 

'1 


The  problem  of  the  steady  cylindrical  source -type  or  sink-type 
flow  has  been  of  interest  to  fluid-dynamicists  for  several  reasons.  First, 
it  is  known  that  the  corresponding  problem  of  an  inviscid  compressible 
fluid  has  an  exact  solution  containing  a limit  line  of  rath^  special  type, 
namely,  the  sonic  circle  (Ref.  1)  . To  the  exterior  of  this  circle  the 
solution  has  two  branch  of  values,  one  has  its  stagnation  point  at  infinity 
(subsonic  branch)  and  the  other  starts  with  maximum  velocity  at  infinity 
(supersonic  branch}.  Both  t the  limit 

line  with  infinite  velocity  gradient.  Therefore  the  viscous  and  heat- 
conductive  effects  arc  expected  to  play  an  important  role  in  continuing 
the  solution  further  inward.  Second,  because  of  its  cylindrical  symmetry, 
this  problem  is  one  of  the  few  nonlinear  flows  in  more  than  one  dimension 
for  which  there  is  only  one  independent  variable,  the  radial  distance. 
Consequently,  the  equations  are  simple  enough  to  allow  a unified  discussion 
of  the  various  effects.  These  are  perhaps  the  reasons  why  this  problem 
has  attracted  the  attention  of  several  authors  (Ref.  2,  3,  4). 

In  the  first  part  of  this  investigation  the  Navier -Stokes  equations 
are  given  for  the  cylindrical  sink  flow  of  a viscous  heat-conducting  perfect 
gas.  The  energy  equation  is  integrated  once  to  give  a first  order 
differential  equation.  Then,  with  some  simplifying  assumptions,  the 
qualitative  properties  of  the  solutions  are  discussed  in  detail  for  the  case 
of  large  Reynolds  number  . The  definition  of  is  , 

where  ^ ^ locates  the  inviscid  sonic  circle  with  sonic  speed  and 

^ are  the  fluid  density  and  viscosity  coefficient  at  /*=^  . These 


basic  properties  of  the  solutions  thus  comprehended  serve  for  a useful 
guide  in  our  final  calculation  of  the  solution. 

In  the  second  part  of  this  paper  the  detailed  calculation  of  the 
solution  is  carried  cut  for  the  case  of  large  ^ . It  is  found  that  there 
is  no  single  expression  available  for  the  solution  uniformly  valid  in  the 
entire  flow  region.  The  calculation  is  then  performed  in  three  different 
regions  characterized  by  the  length  and  the  parameter  ^ For 

f'  > f;  i-  0( ) the  approximate  solution  is  obtained  by  using  the 
PLK  method*.  The  result  fails  to  be  a good  approximation  for  r too 
close  to  . For  -0(Rt.  ^ / 0(&6  , a different 

similarity  rule  for  the  variables  leads  to  a system  of  cylindrical 
transonic  equations  which  governs  the  flow  across  the  sonic  region. 

These  equations  can  be  integrated  analytically  for  each  of  the  different 
order  terms.  The  result  shows  that  the  solutions  belonging  to  the  super- 
sonic branch  all  contain  cylindrical  shock-type  flow  in  this  transonic 
region.  In  other  words,  these  solutions  gradually  deviate  away  from  the 
inviscid  supersonic  branch,  reach  a minimum  near  s ^ and  then 

approach  asymptotically  to  the  viscous  subsonic  branch.  It  is  also  found 

-JU 

that  the  shock  strength  is  of  0(^e  ) while  the  shock-thickness,  of 

results  which  are  quite  different  from  that  of  the  plane  normal 
shock.  Within  this  region,  the  thermodynamic  variables  satisfy  the 
isentropic  relation  up  to  the  order  ) but  deviate  from  it  by  a 

* This  terminology  was  introduced  by  Prof.  H.  S.  Tsien  in  a series  of 
seminars,  held  in  California  Institute  of  Technology  in  1954,  in  v/hich 
the  method  due  to  Poincare^  —•  Lightill  — Kuo  was  diacassed. 


quantity  of  0(^c  . The  approximate  solution  for  r <,  /j  - is 

subsequently  carried  out.  Finally,  the  entropy  variation  of  the  fluid  and 
the  effect  due  to  variation  in  viscosity  coefficients  are  discussed. 

The  corresponding  source  flow  problem  was  previously  solved, 
using  numerical  method,  by  Sakurai  (Ref.  3);  a qualitative  investigation 
on  this  problem  was  later  elaborated  on  many  points  by  Levey  (Ref.  4), 
by  making  use  of  some  conventional  methods  in  nonlinear  differential 
equations.  In  the  latter  work,  the  orders  of  magnitude  of  many  flow 
quantities  of  interest  were  estimated.  The  present  investigation  on  the 
sink  flow  is  not  merely  a special  case  of  the  cylindrical  flow  other  than 
the  source  type,  but  adso  presents  an  improved  method  superior  to  those 
used  in  the  previous  works  (e.  g.  Ref.  2,3,4).  The  powerful  PLK-method 
applied  to  the  outer  region  yields  a set  of  reliable  boundary  values  for 
the  transonic  region  and  thus  enables  all  flov/  quantities  of  interest  to 
be  calculated  quantitatively  in  all  regions.  However,  the  previous  works 
should  be  credited  here  as  to  have  led  this  author  to  a better  understanding 
04  this  problem.  The  author  is  also  indebted  to  Prof.  H.  S.  Tsien  fox 
suggesting  the  problem  and  to  Profs.  M.  S.  Plesset  and  C.  R.  DePrima 
for  their  assistance  on  many  points. 

1.  The  Fundamental  Equations 

Here  we  are  concerned  with  the  two-dimensional  sink  flow  of  a 
viscous,  compressible,  heat-conducting  fluid  ^vith  polar  symmetry.  The 
only  independent  variable  is  the  radial  distance  r from  the  origin.  The 
radial  velocity,  u,  , is  the  only  velocity  component  and  is  always  non- 
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positive  for  sink  flow.  L.et  ^ , f , T , ^ , jn.  , X , ^ denote 

respectively  the  pressure,  density,  absolute  temperature,  coefficients 
of  shear  and  bulk  viscosity,  heat  conductivity,  gas  constant,  specific 
heats  at  constant  volume  and  pressure.  The  momentum  equation  is 


df> 


da 


(1.1) 


and  the  energy  equation  is 


Tr  1';]  ■ 

. 


(1.2) 


The  continuity  equation  can  be  written  in  the  following  form  if  m denotes 
the  sink  strength, 

ZJT  ojir  ~ -m  /I  _ 3) 

The  equation  of  state  is  assumed  to  be  that  of  a perfect  gas, 

^ r . (1.4) 

Equation  (1.  1-1.4)  are  a system  of  nonlinear  differential  equations  for 
four  variables  u , ^ ^ and  T if  ^ ^ , and  are  known 

functions  of  T . 

To  reduce  the  equations  to  nondimensional  form,  the  following 
nondimens ional  quantities  are  introduced: 

f =/•//>,  , f = f/g  . A'A/f-, 


(1.5) 


where  quanities  with  the  subscript  1 are  ficticious  quanties  which  would 
occur  at  the  local  Mach  number  unity  for  nonviscous  and  non-heat- 
conducting gas.  Thus,  with  equal  to  the  ratio  of  specific  heats. 
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assumed  constant  throughout,  the  sonic  speed  a,  at  r-r  is  given  by 


, and 


The  continuity  equation  then  becomes 

^ wr  = / 


(1.6) 


(1.7) 


Here  tv  is  always  positive  for  sink  flow.  The  equation  of  state  is  now 


/.ftf  . 


Using  £q.  (1.7),  and  writing 


the  nondimens ional  form  of  Eq.  (1.1)  becomes 

where  of  denotes  the  inverse  of  the  Reynolds  number 
5f  - A _ 


(1.8) 


(1.9) 


(1.  10) 


(1.11) 


which  will  be  considered  much  smaller  than  unity  throughout  this  paper. 


Stokes  assumption  on  the  viscosity  coefficients  states  that 


X= 


(I.IZ) 


As  this  assumption  does  not  agree  with  observations  for  many  kinds  of 
fluid  (cf.  Ref.  5),  condition  (1.  12)  will  not  be  imposed  on  the  final 
calculation  of  the  flow  field. 

Again  using  £q.  (1.7),  the  nondimens  ional  form  of  Eq.  (1. 2)  can 
be  integrated  once  to  yield: 


(1. 13) 


c 


t 


<_  . 
*7' 


o'  = 

X 


where  (T  denotes  the  Prandtl  number 

(1.  14) 

The  integration  constant  on  the  right  hand  side  of  Eq,  (1.  13)  is  chosen 
as  shown  above  so  that  the  limit  solution  for  vanishing  viscosity  agrees 
at  large  r with  that  of  a nonviscous,  iao-energetic  flow. 

Eliminating  ^ in  Eq.  (1.  10)  by  using  (1.  7)  and  (1.  8),  and 
introducing  j =jlc^F  as  the  independent  variable,  we  obtain: 


(1.15) 


r^/ 


(1.  16) 


The  energy  equation  in  terms  of  j is  now 

Eqs.  (1.  15)  and  (1.  16)  are  the  two  equations  for  two  unknowns  w and  6 . 
The  boundary  conditions  for  them  can  be  determined  by  requiring  that 
they  tend  to  their  respective  inviscid  solutions  as  j~^oo  go  that  these 
two  solutions  can  be  appropriately  compared  later. 

1.  1 The  Nonviscous  Solution 

The  solution  for  sink  flow  of  a compressible  inviscid  gas  can  be 
literally  obtained  by  putting  a - o in  the  above  equations  vrithout 
justifying  the  validity  of  such  a simplification.  Then  the  equations  reduce 
to 

dw 


and 


SI  + ±\  1.(1]  - ±\  = d7 
dj  fvJ 


(1.17) 


(1.18) 


i- 


The  solution  of  this  system  of  equations  is  known  to  be 
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( ^ J 

w <3  y 


f ^ =6/' 


-il.  ZO) 


The  value  can  be  expressed  in  terms  of  stagnation  state  as 

yy-/ 

r.  = - ^ rri  (1.21) 

2/rf,a,  zTta^q,  zita%'<  J 

i Vf> 

where  s _ , and  ^ ^ are  the  isentropic  stagnation  pressure 

and  density.  Equation  (1. 20)  simply  states  the  isoenergetic  and  isei  tropic 
relations. 

This  inviscid  solution  if{r)  given  by  £q.  (1.  19)  is  plotted  in  Fig.  1, 


Fig.  1.  Graph  of  Inviscid  Solution. 


It  gives  no  solution  for  , but  for  , MT  is  a double-valued 

function  of  f On  one  branch  ft  tends  to  zero  so  that  thermodynamic 
variables  tend  to  their  stagnation  values  as  r -►•«>;  on  the  other  branch 
W tends  to  maximum  speed  attainable  \jHi  , and  the  thermodynamic 


variables  tend  to  zero  as  r~»~oo.  They  will  be  designated  as  subsonic 
and  supersonic  branch  respectively.  Both  of  the  branches  terminate  at 
with  sonic  speed  (at  which  the  fluid  speed  equals  the  local  speed 
of  sound).  The  slope  of  the  curve  ^(r)  , 

y 


dr  Yh 


(1.22) 


is  much  smaller  than  unity  for  r»/^  on  both  branches  and  consequently 

viscous  effects  become  comparatively  unimportant  there.  But  as  , 

dtr 

Vf-»>  / , ~ becomes  numerically  unbounded,  and  thus  the  viscous  force 

df* 

near  r = /^  should  play  a role  as  significant  as  those  of  inertia  and  pressure 
forces. 

Now  this  inviscid  solution  will  be  used  as  a guide  to  study  the  sink 
flow  of  a real  fluid  governed  by  Eqs.  (1.  15,  1.  16)  for  large  '.values  of  He  , 
in  the  sense  that  it  is  assumed  that  the  limit  of  the  viscous  solution  for 
vanishing  viscosity  approaches  the  inviscid  solution  as  r-9-o*>  , for  both 
subsonic  and  supersonic  branches.  By  continuing  these  viscous  solutions 
backward  in  r where  viscous  effects  become  more  and  more  prominant, 
it  is  expected  that  the  real  fluid,  affected  by  viscosity  and  heat  conduction, 
will  flow  across  this  fictitious  sonic  circle,  which  is  a limit  Une  when 
viscosity  is  neglected. 

It  may  be  remarked  here  that  the  equations  foz^  source  flow  of 
real  fluid  can  be  obtained  from  (1.  15)  and  (1.  16)  by  changing  the  sign  of 
terms  with  factor  d(  if  W again  represetzts  the  absolute  value  of  the 
radial  speed,  normalized  relative  to  (cf.  Ref.  4).  Hence  the 

inviscid  solutions  for  source  and  sink  flow  are  identical,  but  their  respective 


viacouB  solutions  will  be  shown  later  to  have  quite  different  features  for 
ail  r . 


2.  Properties  of  the  Solution  Curves 


2.  1 Approximate  Differential  Equation  in  the  "Phase  Space". 

In  order  to  study  the  qualitative  properties  of  the  solution  curves, 
several  assumptions  will  be  introduced  in  this  section  to  simplify  the 
analysis  while  most  of  the  important  features  of  the  original  system  will 
still  be  maintained.  The  Prandtl  number,  O'  , is  assumed  constant 
because  yi  and  'X  have  almost  the  same  dependence  on  temperature. 
In  this  section,  /t  is  also  taken  to  be  constant  so  that  ' • When 

the  complete  solution  is  calculated  later,  however,  the  asstunptions 
introduced  here  become  unnecessary. 

Equation  (1.  13)  can  be  integrated  when  the  Prandtl  number 


(T  = 


H-4-^  i 
Z(ti-l) 


(2.1) 


(under  Stokes  assumption,  4=  » tken  0"=  ^ ~ ^ 

integral  is 


(f  ^ 


(2.  2) 


where  A is  the  integration  coi^stant.  The  value  chosen  above  for  (f  is 
actually  not  far  from  experimental  data  { ^ =.  o.72  for  air  at  standard 

</e 


condition).  As  O'  only  appears  in  the  coefficient  of  the  derivative 


in  Eq.  (I.  lb),  it  follows  from  the  theory  of  differential  equation*  (Ref.  6 
p.  142)  that  the  solutions  and  all  their  derivatives  will  be  continuous  in 
O'  for  vf  > O , _oo<^<oo  . Thus  the  assumption  of  choosing  this 


'1  f*  T-iv- 


M 


i 

>■* 


V. 


particular  value  ox  ff'  would  merely  lead  to  simplification  of  analysis 
rather  than  material  change  of  the  solutions.  If  we  further  require  by 
physical  argument  that  the  deviation  from  the  iso-energetic  relation 
expressed  by  the  term  with  the  arbitrary  constant  A v/ill  not  overwhelm 
the  right  hand  side  terms  for  , we  may  assume  that  A = o . This 

restriction,  however,  can  again  be  relaxed  when  the  complete  solution 
is  discussed  later.  It  will  then  be  shown  that  A is  indeed  of  the  order 
o(e()  , Thus  the  particular  solution  with  A=o  would  still  provide  a good 
approximation  to  the  complete  solution. 

With  A=0  , the  energy  relation  becomes 

9 = - tzL  (/i-4.6(^)vr^  . (2.3) 


Introducing  this  equation  into  £q.  (1.  15)  and  eliminating  the 
explicit  dependence  on  'rj  by  the  substitution 

U-_ 

~ di^  ~ dr  * 


(2.4) 


we  obtain 


avf*y  ^ / ^[/  - - iV  t -f-  ■ - O (2.  5) 


where 


} f = 


a - 


(2.5a) 


tH'-  ' r - uf  ' ~~77I^T 

The  variable  V is  closely  related  to  the  fluid  velocity  gradient.  Since 
the  terms  acx  and  (a-i)fA  in  the  brackets  are  merely  corrections  to 
constant  coefficients  of  O(l)  , the  properties  of  Eq.  (2.  5)  would  not  be 
altered  if  we  had  neglected  these  terms  in  order  to  simplify  further 
algebra.  Thus  the  approximate  differential  equation 

~ = o (z.fi) 

jittr 
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ia  the  phaee  space  (*>'  ,V)  is  expected  to  exhibit  all  important  features 
of  ^he  original  system  Eqs.  (1.  15  - 16)  for  /*  = >'  . The  equation 
similar  to  (2.  6)  wac  derived  by  Sakurai  (Ref.  3)  and  later  was  discussed 
in  detail  by  Levey  (Ref.  4)  for  source  flow  in  a real  fluid. 

2.  2 Properties  of  the  Solution  Curves  in  Phase  Space  . 

Equation  (2.  6)  is  nonlinear  and  cannot  be  integrated.  However, 
several  important  features  of  the  solutions  can  be  readily  seen  by  study- 
ing the  properties  of  the  vector  field  defined  by  Eq.  (2.6),  such 

as  the  type  of  its  singular  points,  the  curves  of  zero  slope  and  zero 
curvature  together  with  some  obvious  isoclines. 

(a)  The  curve  of  zero  slope;  the  inviscid  solution. 

Let  be  the  curve  on  which  ~ given  in  Eq.  (2.  6)  vanishes, 
C,  is  then  given  by 


(2.7) 


which  is  also  the  inviscid  solution  in  w-V  plane.  The  function 

given  by  Eq.  (2.7)  has  a simple  pole  at  vv'a/  (the  fictitious  sonic  circle). 


and  two  zeros  at  o and  W = ^ 


which  correspond  respectively 


to  the  subsonic  and  supersonic  branch  at  r—oo  . Near  the  origin,  V(w) 


has  the  following  power  series  expansion 


i- (/-^)w  ‘ (2.8) 

; 

which  starts  from  W = 0 with  slope  unity.  Near  the  point  yr  = ^ , 

the  expansion  of  Y(^)  is 
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where 


= w - ^ . 

/ 

This  branch  of  the  curve  starts  from  ,o)  with  the  slope  — 

t 

The  curve  divides  the  infinite  strip  ^ ^ into  regions  of 

positive  and  negative  slope  as  shown  in  Fig.  2. 


Fig.  2.  Regions  of  Positive  L.nd  Negative  Slope;  Isoclines  . 


In  this  strip,  ^ 


given  in  Eq.  <2.  6)  becomes  infinite  only  on  {i)  V'so, 

« ' 

o<w<j5‘*  , and  (ii)  w=o  , Y^O  . Besides,  iV  = ^ ^ , Y^o  is  also 

an  isocline  on  which 

<JV 


dl/f  PT 

(b)  Properties  of  the  singular  points. 

The  only  singular  points  of  Eq.  (2.  6)  are  and  (/^^/ 

The  origin  is  a singular  point  of  higher  order.  But  if  one  sketches  the 
vector  field  defined  by 


^ = w(/-fvr)  ~V  ( /~nr) 


rigin  is  a saddle 
One  of  diese  is 
has  no  physical 


along  a simple  closed  curve  dd)  in  the  neighboSWood  of  the-origin  with 
the  origin  in  its  interior,  one  finds  that  the  PoinCv^^j^;'  i,idex  (c£.  Ref  7 
p.  45)  of  this  singularity  is  equal  to  - / . Thus  the\' 
point  thi'ough  which  only  two  solution  curves  may  past 
v:  - o 'vhich  either  yields  a trivial  solution  (V^^o)  o' 

meaning  { y % o)  . The  other  solution  curve  starts  from  origin  with 

slope  equal  to  V-  / (which  coincides  there  with  the  inviscid  -^KVolution)  and 
thus  represents  the  only  possible  radial  sink  flow  with  stagna'i^^\on  r*oo. 
By  substitution  of  a power  series  into  £q.  (2.  6)  (or  by  ordinart^^^ 
the  asympi^otic  value,  for  small  o(  , of  this  solution  near  the  o; 
found  to  be  " 

V ^ w [(/-/!)- ajv/  y-  [(/-/>) -« 


iteration), 
igin  is 


y-  [('if)  -0i(U  - i6^  : }f)t/0(K  y-  • 


^\(2>  1,1) 


The  .point  ,o)  is  » singular  point  of  regular  type.  In  the 

J 

neighborhood  of  j;3W-^^=o,  V«o  , Eq.  (2.  6)  becomes 

dV  _ j2.  il) 

otV  y- 


</e 

b * 

in  which  P and  ^ anish  like  x y-  V 
singularity  (Ref,  7 pp.  37-44)  is 
(i)  a nodal  point  if  oi  /■ 


('-**■  y 


ae  X,  ]f-^o  , Thus  the 
, (for  air,  ^ ^ ^ 


c-P) 

(ii)  a spiral  point  if  or  > ~7rj~ 

ep 


(2.  12) 


As  the  problem  will  be  confined  to  the  case  o(  « f , we  sh2dl  only 
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consider  this  singularity  to  be  a nodal  point,  (which  changes  to  a saddle 
point  for  source  flow,  cf.  Ref.  4 ).  All  solution  curves  passing 
through  this  point  will  have  at  this  point  two  distinct  slopes  which  can  be 
calculated,  from  the  secular  equation  of  Eq.  (2. 11),  namely. 


■X. 


-of 


2^ 


O 


(2.13) 


which  has  two  unequal  positive  roots 

^ + ^ - ^ ■ <2-  “»> 


From  these  two  eigen-values  two  eigen-vectors  associated  with  Eq.  (2. 11) 
at  (z^o , V^o)  can  be  obtained.  By  using  this  result  it  can  be  shown 
that  the  solution  curves  passing  through  (x  ^ V"0)  have  near  this 
point  the  following  parametric  representation 


(2i  14a) 


with  the  slope 

(fy-  _ 

efx  • t't/Jt 


(2. 14b) 


Since  >o  , V$X-*o  as  /->•  -®<»  . It  also  follows  from  the  above 

equations  that  there  are  infinite  number  of  solution  curves,  corresponding 
t!i  arbitrary  C,  and  which  have  the  asymptotic  value 

near  (z^O,  y=^0)  ; (2.15a) 
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Vi-*  * 


and,  in  addition,  there  is  another  solution  curve  passing  through  this 


point  , of  the  value 

V = i—  - — ) 

S(  t-fi  ' 


(x=c  o,  V = o)  . 


(2.  15b) 


The  first  group  of  solutions,  given  by  E^.  (2.  15a),  have  the  same  limiting 
value  at  rsoo  as  the  inviscid  solution  and  hence  represent  the  many 
possible  sink  flows  starting  with  maximum  velocity  w - ^ ^ at  r=  oo  , 
while  the  solution  given  by  Eq.  (2.  15b)  is  physically  irrelevant.  The 
asymptotic  value,  for  small  or  , of  this  physically  significant  solution 
near  the  point  is 

2 1 


O'fif/ztp 

y.  C iit 


(nea.T  z ^ * so  ^ 


(2. 16) 


where  C is  an  arbitrary  constant.  The  last  term  follows  from  Eq.  (2.  15a). 
(c).  The  curve  of  zero  curvature. 

The  second  derivative  of  the  solution  V can  be  deduced  from 


Eq.  (2.6), 


3 3 

<X\/fY 


= zV^-w(/^fivr)V  -f- ^ (f-vr)-if(f-/}vr)^  . 


(2.17) 


Let  be  the  curve  on  which  vanishes,  the  equation  for 

is  then,  except  where  vfV  is  zero,  given  by; 

2V^-vf(/^^W^)V  -f-  ^ =o  (2.18) 

The  function  satisfying  this  equation  has  the  following  properties: 

(i)  has  only  one  real  value  for  eithex 


'/3 


04  or  ; 


(2.  19a) 


lA  has  three  real  values  for 


y.  has  two  equal  real  values  at 

M 


{2.  19b) 


•'■f^  ■ (2.19c) 


(ii)  The  curve  starts  from  the  origin  with  slope  equal  to 

unity  and  has,  in  the  neighborhood  of  the  origin,  the 

following  expansion: 

j z ^ 

w -f-  [(/-^) fj3)t-fo(  J W 

^ \('-f)-ofO-f  -/^/3  . { 2.  20) 


(iii)  The  curve  crosses  w*/  at 

y^O)  ~ (i^) 0(<x^  )'^ 

with  the  slope 


^ = 1 0-f)  r / 1 


(Iv)  When  w"  = y?  ^ ^ - <7 


Thus 


/^-^o  and  \=f^ 

?v  double  point  at  oj  , where  the  curve  has  two 

different  slopes  [i  + 0(c()J  and  ^ 0(b()J  . Near 

this  point  these  two  branches  of  the  curve  has  the 

following  expansions: 


(2.  21a) 

{2.  21b) 
has 


W 


1 


■■■  ] 


(2.  22a) 
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(2.  22b) 


(v)  The  slope  of  the  curve 

4Y^-3W'(/^/tM-^)V  V-  ^ (f-Yf*)(>~fin^^) 


(2.23) 


becomes  infinite  at 

(vi)  For  «•»/.  j fS  . (2.24) 

The  curve  divides  the  (*r~Y)  plane  ?nto  regions  of  positive 
and  negative  curvatures  as  shovm  in  Fig.  3. 


Fig.  3.  Regions  of  Positive  and  Negative  Curvature. 
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1- 


/ 

i 


(d).  Sketch  of  solution  cufves  Vv'-.'')  ; definition  of  cylindrical  shock. 
The  above  discussions  on  isoclines,  types  of  singularities,  regions 
of  positive  and  negative  slope  and  curvature,  enable  the  solution  curves 
to  be  sketched. 

First  let  us  consider  the  solution  curve  starting  from  the  origin. 
Comparison  of  Eqs.  (2.8),  (2.10),  and  (2.  20)  shows  that,  for  oc  and  w 
small, 

V (^)  > V(w)  > V (w)  , (2.25) 

the  difference  (V^  -V)-sjO(c(W^)  while  (V-\^  )~  . Thus  the  solution 

curve  V(*f)  lies  in.  between  ^ and  where  the  slope  and  curvature 
of  V(w)  are  both  positive.  A careful  study  of  the  slope  of  VW  and 
(cf.  Eqs.  2.  6 and  2.  23)  indicates  that  the  solution  curve  lies  above  for 
w>o  . Hence  V(<^)  is  a monotonic  increasing  function  of  w with 
increasing  slope,  passing  through  w = / between  the  points  V = and 

V = , and  finally  ending  up  at  w = with  the  slope  ^ =.  — £ 

' aiir 

(cf.  Fig.  4). 

As  previously  shown  in  Eq.  (2.  16),  there  are  infinite  ntunber  of 

• — - 

solution  curves  starting  from  w = y?  ^ , V=o  with  the  same  slope  ■■  ^r.  • 

.X 

However,  for  and  o(  both  small  enough,  comparison  of 

Eqs.  (2.9),  (2.  16)  and  (2.  22a)  again  shows  that  all  these  solutions  satisfy 

^ / 

for  small  negative  , the  following  inequality 


V(vr)  < V.  in)  < V (w) 


(2.  26) 


as  shown  in  Fig.  4. 

/ 

As 

decreases  with  increasing  slope  until  it  intersepts  ^ with  a positive 


• 18. 


Fig.  4.  Sketch  of  the  Solution  Curves  in  the  Phase  Space  . 


slope.  For  further  decrease  in  w , the  curve  V should  lie  above  a 
straight  line  with  this  positive  slope  at  the  point  of  interception  because 

</V 

is  positive  in  this  region.  Hence  the  solution  curves  will 
eventually  meet  (*  with  zero  slope,  and  from  thve  on,  for  further 

decrease  in  , V increases  from  negative  values  and  later  crosses 

/ 

y s o with  infinite  slope  at  some  point  in  between  = o and  vV=  y?  ^ 
as  can  be  shewn  by  the  method  of  bounding  curves  and  will  be  made 


r,/-- 


m- 

\ 


explicit  in  our  later  calculation.  Further  extension  of  these  solution 
curves  shows  that  V incr?ases  with  increasing  w and  finally  approaches 
asymptotically  to  the  subsonic  branch  solution  which  starts  from  the  origin. 
There  is  a particular  value  of  the  integration  constant  (in  Eq.  2.  16),  say, 
c = q,  < o , for  which  the  solution  curve  finally  ends  up  at  the  origin  with 
infinite  slope.  For  c < , the  solution  ceases  to  have  physical  meaning. 

On  the  other  hand,  the  solution  curves  for  c>  have  a very  interesting 
feature  that  these  viscous  solutions  all  exhibit  the  transition  process  from 
the  "inviscid  supersonic  branch"  toward  the  "viscous  subsonic  branch". 

L.et  us  consider,  in  particular,  the  solution  with  c - o . It  first  inter- 
cepts at  , BA',  and  then  crosses  V=o  at  . Since 

>/  and  ^ (S’S  will  be  shown  later),  the  flow  between  these  two 
states  may  thus  be  defined  as  that  of  a "cylindrical  shock".  * Inspired  by 
the  result  obtained  in  Eqs.  (2.21),  we  see  that  the  equation  governing 
such  a cylindrical  shock  flow  can  be  approximated  by  the  following 


similarity  transformation 

n - / i-  « w 


(2.27) 


J 

]\ 


This  terminology  is  adopted  by  both  Sakurai^  ^ and  Eevey^  ^ to  describe 
such  type  of  flow.  The  term  "shock"  is  borrowed  from  its  conventional 
meaning  to  indicate  the  transition  from  one  branch  to  the  other,  though 
the  transition  is  rather  different  from  that  occurring  in  a plane  normal 
shock.  Perhaps  this  terminology  relates  closer  to  the  conventional 
meaning  of  a shock  for  the  constant  C slightly  greater  than  c„ 

(cf.  Fig.  4),  because  then  the  jump  in  w and  ^e  slope  ^ in  transition 

become  greater  and  the  position  of  transition  is  farther  out  from  r = /; 

(cf.  Fig.  5).  But  since  there  is  no  adequate  criterion  to  distinguish 
one  from  another  value  of  C , we  shall  retain  this  name.  Another 


terminology,  the  dissipation  layer,  is  suggested  by  Prof.  H.  S.  Tsien 
to  void  this  ambiguity  and,  in  addition,  to  stress  thb  importance  of 
viscous  effects  in  this  layer. 


which  will  render  all  terms  in  £q.  (2.  6)  equally  important  in  this  flow 


region,  that  is,  for  vanishing  a , 

V = ZyfV  (/-p) 

dff 


(2.  28) 


Since  this  equation  governs  the  flow  of  both  branches  near  the  sonic 
speed,  it  may  be  called  the  "equation  for  cylindrical  transonic  flow", 
(e)  Sketch  of  solution  curve.<3  in  physical  space. 

From  the  definition  of  V , - -rr  > we  have 


(2.29) 


where  the  integral  stands  for  an  indefinite  integral  and 

SO  that  fj  tends  to  its  inviscid  solution  for  j large. 

Now  for  the  subsonic  branch,  , hence  y is  a 

monotonically  decreasing  function  oi  ■ w . Moreover,  for  same  value  of 
w , V(w)  is  less  than  its  corresponding  inviscid  solution  (cf.  £q.  2,  25). 


(2.29a) 


Hence,  from  Eq.  (2.2  9), 


. (*v)  < 


-(2.30) 


In  other  words,  at  every  y , is  slowed  down  from  its  inviscid 


value  due  to  the  viscous  effect. 


For  the  supersonic  branch  starting  from  w =f  ^ . V(^)  <o 

I 

for  /T  < W $ ^ , hence  in  this  interva?  y is  a monotonically  increas- 
ing function  of  W . At  yf-yfg  » ^ ®ayj. 


(^)  =_V(lAr)=0 


and,  from  Eq.  (2.  6) 


Wht4,o  ^ 


i 


t 


therefore  W(^)  = is  the  only  minimum  of  w on  this  branch.  For 

J ^ ^ decreases  vrith  increasing  ^ , Furthermore,  because  of 

the  relation  given  in  £q.  (2.  26),  we  have 

7 for  , V^o  . (2.31) 

A«  /«y« 

The  above  properties  of  the  solution  enable  the  solution  curves  to 
be  sketched,  as  shown  in  Fig.  5.  As  the  constant  c decreases  from 
zero,  the  minimiun  value  decreases  and  the  jump  in  W increases  while 
the  jximp  takes  place  farther  upstream. 


c 


1 . 

*>■ 
r'  ■ 


Fig.  5.  Sketch  of  the  Solution  Curves  in  the  Physical  Space. 


c 


It  follows  from  Eq.  (2.  27)  that  if  we  introduce 

'/s^ 

j ^ cK  ^ , >^=/v‘crw. 


then  Eq.  (2.  28)  becomes 

A 
¥ 


<^0 


= -Z'tf  ^ -hi /-f) 


which  is  the  equation  governing  the  flow  near  w s / , j 


= 0 
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- (2.32a) 

_ (2.  32b) 
Integrating 


(2.  33) 


this  equation  once,  we  obtain 


y.  ^ = (/~f)^  f 


Consi. 


This  equation  will  be  integrated,  and  discussed  in  detail  in  cur  final 
calculation. 


3.  Calculation  of  the  Solutions  by  Using  PUC-Method* 

In  this  section  we  shall  calculate  Hfl'tj)  , 6(j)  governed  by  the 
original  system  of  Eqs.  (1.  15  - 1.  16).  Throughout  this  section  jH  will 
again  be  assumed  constant  so  that  /L-t  , but  no  restriction  will  be 
imposed  on  o'  and  ^ . Consequently  Eqs.  (1,15,  16)  become 


where 


o(  = ^ 


y 77T 


Using  the  PEK-method,  as  described  below,  the  generalization  to  the 
case  yC  -^Ct)  presents  no  particular  difficulty  (cf.  J 7). 

It  was  seen  before  that  even  for  a simplified  version  of  these 
equations,  such  as  Eq.  (2.  6),  conventional  perturbation  method  merely 

i 

leads  to  asymptotic  solutions  for  small  or  near  <k-0  or  w'  ^ f ■*  (cf. 

Eqs.  2.  10,  2.  16)  because  coefficient  of  w'^  does  not  diminish  as  /?-«■<>»  , 

This  asymptotic  resvilt  fails  to  be  a good  approximation  to  the  required 

$ 

solution  as  W~  deviates  farther  from  W=o  or  and  becomes 


* Cf  the  footnote  in  the  Introduction. 
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i 


almost  useless  for  calcolatioas  near  vr  -/  < Now  let  us  resort  to  the 


PLK -method,  which  is  in  essence  to  expand  the  solution  in  terms  of 
power  series  in  or  with  coefficients  as  undetermined  functions  of  a 
parameter  ^ , 

0)  2 a) 

^ -h  o(  v^(^)  -f- a w (^)i-  • ■ ■ , 

fa>  0>  2 (2) 

e(^)  ^Q($)'ho(6(^)i-o(  G (^)f-  . ....  .. 


(3.4) 


The  need  of  a parameter  ^ to  represent  the  soluti'^n  and  that  W starts 
with  the  term  ^ are  clearly  suggested  by  our  previous  discussions. 
Substituting  these  expansions  into  Eqs.  (3.  1)  and  (3.  2),  noting  that 
ditf  _ wr  d6  _ S 

HL  r 

^ f ' ^ ‘ 7'  ' df  ' (y'f 

where  prime  stands  for  ~ , and  then  equating  equal  powers  of  or  , 


jf. 

Y > 


(3.5) 


we  obtain  the  zeroth  order  equations  as  follows: 


(3.6a) 


(3.6b) 


If  we  choose  different  from  zero,  then  we  have  the  zeroth  order 


solution 


As)  = ^ ¥ I 


(3.7a) 


(3.7b) 


The  integration  constant  in  Eq.  (3.7b)  has  been  so  chosen  that  when 

(0) 

, 0 - vf  and  J (vf)  agrees  with  inviscid  solution. 

The  coefficient  of  in  the  expanded  equations  gives  the 


following  first  oxder  equations 

d‘’i-  (Y->) 0 vr'^ --  - ~ + ~f. 


r 


(3.  8a; 


where 


(3.  8b) 


^ , /=34«y,,  ,3.8c) 


-i  H-i  ' / ( tn  ’ jr  *-• 

By  substituting  Eqs.  (3.  7)  and  (3.  8a)  into  (3.  8b),  the  terms  in  the  curly 
bracket  can  be  rearranged  to  talce  the  follov/ing  form 

If'  ^ m ^ 


\bB~e  \ 

L r ■»  ^ 


and  finally  Eq.  (3.  8b)  turns  out  to  be 


r (0)^0)  0 

« r - « 

^ ^ f h 

~CL)0-fi) 

i 

L/  / 

J L 

• Y 

- r /^_6) 

1 

•f 

X— 

» 

Vr 

> 

(3.  9) 


r y I 


-ji  % 


C 


Now  in  order  to  ascribe  to  ' a possibly  lowest  order  singularity  at 
^ = / to  improve  the  convergence  of  the  series,  we  decompose  Eq.  (3.  9) 
as  follows 


^ L (1-^^)^  -I 


(3.  10a) 


(3.  10b) 


The  solution  of  Eqs.  (3.  10a, b)  is  thus 


7"w  - f'-o[7fr«  , 


{3.  11a) 


(3.  lib) 


where 


0-a-)0-p) 

A 2 , 

¥ 


(3. 11c) 

In  the  above  first  order  solution.  and  have  singularities  at  ^ ~i 

O' 

of  the  same  order.  After  having  determined  vy  ' , 6'  is  then  given 
by  Eq.  (3.  8a). 

Proceeding  in  a similar  manner  to  obtain  the  second-order 

z 

equations  by  equating  terms  with  of  , we  find  that  the  resulting  equations 

(21 

possess  the  solution  of  quite  lengthy  expression,  in  which  y (0)  starts 

with  the  term  folloy/ed  by  terms  of  o(-, — ) while 

(2)  ^ \ 

W (^)  still  can  be  made  to  be  of  C/(  j . Therefore  the  final 

solution  can  be  expressed  parametrically  as  follows: 


'^1 

n 


-Zb- 


#>' 


r 


-,3  -rz'.f‘\  fO(^)  (3,13.) 


.it 


i: 


^ at .±Ii-/^Jl-Jl  ^ o( 


■h  c 


r I r^/ 
2 


(-fit 


0-^^f  d.f 


_^L_) 

,-^7) 


where  A , ^ ^ are  constants  given  in  Eqs.  (3.  He),  (3.8c)  and 


C *0  for  the  solution  starting  from  k = o , 


C is  arbitrary  for  the  solution  starting  from  ^ » 
The  value  of  6 is  given  by 


/ 

r 


^ «(Y-o[a  ^ 3 ^ ^ f'  ( 


~T^  ^ ^ “7->^  ■ 


(3.  12b) 


(3. 12c) 


(3.  12d) 


( 


Several  interesting  features  of  the  above  solution  may  be 


mentioned  here,  (i)  $ = o for  <T  = 


2 Of/) 


( = 


if  = cf.  Eq.  3.8c). 


With  the  value  of  <r'  and  4 lying  in  the  experimental  range,  c is  still 
a amall  number.  Consequently,  the  variation  in  ^ and  4 only 
contribute  a small  correction  to  the  coefficients  of  tP(7J  in  the  solution. 
This  fact  confirms  our  previous  statement  in  § 2.  1.  (ii)  By  substituting 
£q.  (3.  12)  into  Eq.  (2.  2),  it  can  be  found  that  the  arbitrary  constant  A 
is  of  0(°(^)  . (;,ii)  The  most  important  property  of  the  above  solution 
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■A 

\Kr. 


I 


^11 


( 


is  that  they  do  net  provide  an  approximate  solution  with  uniform 
accuracy  over  the  interval  of  ^ such  that  o i.  vr  ^ . As  ^ 

approaches  unity,  the  higher  order  terms,  especially  in  , become 


more  important  relative  to  the  zeroth  order  term.  More  precisely,  the 


'/3  73  ~2 

solution  is  good  only  for  0^  f ^ !—K<^  and  H- K<x  4^  ^ 


'/j 


K 


Ys 


being  a constant  of  0(>)  . At  ^ ~ i Ko(  , all  terms  in  the  expression 


for  y(^)  become  of  the  same  order,  ) ■,  but  the  convergence  can  be 

made  sufficiently  rapid  by  an  appropriate  choice  of  the  value  K 

In  the  subsequent  calculation  of  the  solution  through  the  transonic 
region,  we  shall  only  consider  a particular  solution  with  C = o , c^.o 


in  £qs.  (3.  12).  Furthermore,  it  has  been  found  convenient  to  take 


K=2(n0  . With  this  valueof  K , we  obtain,  from  Eqs.  (3.  12a,b), 

the  following  result: 


(i)  supersonic  branch, 
at 


(ii)  subsonic  branch 


^ s 2.28/ (Yh)  or  and 


Jw  rt'* 


(3.  13a) 


at  , 


^ 3/j 

y = /./66  (Yh)  o(  and 


.(3  13b) 


if-hO' 


These  values  will  serve  for  the  boundary  conditions  imposed  on  the 
transonic  solution  to  be  obtained  below.  That  the  PJLK  method  is 
powerful  in  solving  this  problem  can  still  further  be  stressed  by  the 
following  argument.  As  the  first  order  term  in  the  expression  for  W(^) 
(cf.  3. 12a)  is  quite  unimportant  in  the  aforementioned  regions  for  0 , 
one  perhaps  would  try,  instead  of  £q.  (3.4),  a simpler  expansion 


(O) 


0) 


j(w)  = y (f/J  + c{  tj  (w)  f • • 


(3.  14) 
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J' 

'41 

.3 1 


i^nd  a similar  expansion  for  ^ in  terms  of  ^ . It  can  be  shown  th?t 

(o) 

the  above  expiansion  will  yield  a solution  in  which  ^ is  identical  to 

0) 

inviscid  solution,  but  y (<^)  has,  in  ?.ddition  to  a simple  pole  at  W=/  , 

' ^ $ 

a pole  and  a logarithmic  singularity  at  Vf  - f . Consequently  the 
assumed  expansion  (3.  14)  becomes  invalid  for  r large  on  the  supersonic 
branch,  and  thus  leads  to  an  erroneous  result. 


4.  The  Solution  in  the  Transonic  Region 

To  obtain  an  approximate  solution  in  the  transonic  region,  as 
discussed  in  § 2.  2(d),  (e)  (cf.  JElqs.  2.  27,  28  and  32)  and  also  as 
guided  by  the  boundary  conditions  Eq.  (3.  13),  first  we  distort  the 


independent  variable  by  the  transformation 

7 = « ^ , 

and  then  expand  vr,  9 into  the  form 

= / 7<-  t f • • • 

e(^)  = / -h i- o('^ i^or9%)  /<•••• 

With^.o  (~^<r(>/;^)=/),  Eqs.  (3.  1)  and  (3.  2)  reduce  to 

I / - vy  1 ^ y.  w / / - (a-Oof]  y/  ] = or  W* 

L J J 

e = — - Iti  (n-ho()vr  -h  0(o(‘) 


(4.  1) 


(4.  2) 


(4.  3) 


(4.  4) 


(4.5) 


Substituting  Eqs.  (4.  1)  - (4.  3)  into  Eqs.  (4.  4)  and  (4.  5),  we  obtain  the 


first  order  eouation: 


(4.6) 
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and  the  second  order  equation: 


J / /.«  /^»  / 3 (f\ 

^ y.  2 = JLfy^  ) . 

‘Js  di  ■ ' 


(4.7) 


The  correction  due  to  the  terms  with  constants  a and  ^ enters  only 

(3)  J3) 

in  w , o and  higher  order  terms. 


Now  Eq.  (4.  6)  can  be  integrated  once  to  yield: 


<0  z 

^ + \tf  = ~£_  % -h  D, 

d§  Hi-!  ^ 


(4.8) 


usiyig 


v/here  the  ccnotant  £)  can  be  determined  by^the  boundary  condition  (3.  13); 
the  result  is, 


= O.OS-3  (ti-l) 

B,  =-0.ct  (t-H) 


-2/3 


■2/3 


for  the  subsonic  branch, 
for  the  supersonic  branch. 


(4.9) 


It  is  convenient  to  rewrite  Eq.  (4.  8)  in  the  following  form 


dx,  <r 


where 


Now  Eq.  (4.  10)  is  of  the  Riccati  type,  which,  by  the  transformation 

,u>  = ± ^ , 

</  vdx, 

can  be  reduced  to  a second  order  linear  equation: 

dtr  . 

The  solution  of  this  equation  for  (xtZ,)yo  is 

'V  ~ M I -f- , £=  ± fxy-JTy/*'  , (4.13) 

where  I(x)  is  the  modified  Bessel  function  of  the  first  kind,  and  M , /V 


(4. 10) 


(4.  10a) 


(4. 11) 


(4.  12) 
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are  the  integration  constants.  By  using  i:jq.  (4.  11),  the  solution  oi 
Eq.  (4.  10)  is  then 


/ , , It,  (z)  -f-  N 

y(z)  ={-^/ 

^ I^A2)+NI,(2) 

-73  A 


where 

£ = A (zfx.,)  , and 

The  constant  /V  can  be  determined  by  using  the  condition  (3.  13). 
The  continuation  of  £q.  (4.  14)  into  the  region  x&z,  <<?  is 


provided  by 


i/j  '/j 


z/*  •Is 

2 ^ 


•/s  '/3 

2 I,  U)  = ^ Jj,  (0 

~3  ~S 

Z (Z)  - ^ 'i2/S%^ 


wnere 


-3ltL/z  J. 


(4.  14) 


(4. 14a) 


(4. 15) 


Consequently  Eq.  (4,  14)  becomes 

pff)  = (^if{  *'  i < o . (4  14, 

To  discuss  the  above  solution,  we  first  note  that  the  inviscid  solution  in 
this  transonic  region  is 

y = X (4.  17) 

which  has  two  branches  for  x >0  and  gives  no  solution  for  x < o . 

Now  before  we  determine  the  value  of  ^/  for  the  corresponding  viscous 
solutions,  we  may  also  note  that  the  genezed  solution,  given  by  (4.  14) 
and  (4.  16),  is  a semi-trancendental  function  of  and  the  second 

integration  constant  //  . It  can  be  shown,  from  the  properties  of 
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at  large  Z , that  in  i^q.  (4.  14) 


/V  = -/ 


and  ^ has  a simple  pole  at  a certain  finite  z for  /V<'-/  (which  is 
of  no  physical  significance).  This  result  shows  that  the  viscous  solutions 
tend  to  their  respective  inviscid  values  for  z large  in  a manner  which 

t 

implies  again  that  iB  a.  nodal  point  (admitting  more  than 

one  value  of  <V  ) while  (vf=o,  y=<>^)  is  a saddle. point  (admitting  only  one 
value  of  N J.  However,  for  xi-z^  <o  , Eq.  (4.  16)  shows  obviously 
that  bas  infinite  number  of  isolated  simple  poles  at  where 

the  denominator  vanishes.  Since  the  properties  of  the  solution  curves 
in  the  (tf  ,V)  phase  space  exhibit  no  such  singularities,  the  solution  (4.  16), 
therefore,  presents  a good  approximation  to  ftie  read  flow  only  for  ^ 
lying  in  the  interval  0(.^  where  is  the  first  pole  and  S is 

a positive  number,  appropriately  chosen  such  that  is  not  yet  too 

large  to  void  our  approximation  (4.  2). 

(0.  / 3 

Having  obtained  the  first  order  solution  W =(-£-)  ^ , the 

second  order  equation  (4.7)  can  be  then  integrated  to  yield 

(2)  , 

w (§)  - e ^ (4.18) 

where 

and 

(z)  <'K 

It  is  obvious  that  w (f)  is  bounded  wherever  vf  (^)  is  bounded. 
Consequontly,  the  approximation  is  good  even  if  we  only  take  the  first 
two  terms  in  (4.  2)  and  (4.  3). 
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In  order  to  obtain  some  numerical  results,  we  first  determine 
the  value  N in  Eqs.  (4.  14),  (4,  16)  by  using  conditions  (3.  13)  to  obtain 


N-  -0.S8S 


N - -t 


for  the  supersonic  branch, 
for  the  subsonic  branch. 


(4.  19) 


With  these  values  of  D,  and  N icf.  4.  9 ar-’i  4.  19),  the  solutions  are 
plotted  in  Fig.  6 (by  using  tables.  Ref.  8)  from  which  several  interesting 
results  can  be  deduced  as  follows; 

(i)  For  the  supersonic  branch,  the  transonic  solution  starts  from 
point  A (cf.  Fig.  6)  with  the  coordinates 


(ii) 


= ! i-  f.6 


/j  >/3 
a( 


^/s 


i z.9(Uiy<x 


(4.  20) 


y 

After  the  solution  curve  passes  through  a point  of  inflection  G 
and  then  crosses  the  line  ^~o  at  point  £ (cf.  also  Fig.  4) 

where  =.  i.oz  Vg  ^ ^ 

{—]  =1  ~0’Ss{~)  cK  ) , it  reaches  a minimum  when  it  inter - 

cepts  the  curve  ^ = .c  at  C where 

/ -2  '/j 


(4.21) 


It  then  increases  from - VT  = to  w = / at  point  D where 

That  is,  there  is  an  expansion  wave  following  the  cylindrical  shock. 
The  thickness  of  the  cy’indrical  shock,  as  defined  in  § Z,  2(d),  i;> 


Ay-=.y^-y^  =■  Z-7  o<  (n  /*-  space  , Ar  Ay^  (4.22) 


across  which  the  velocity  has  a "jump" 

jy-3  7j 


AW  = W-W,  «' 


(4.  23) 
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Combining  £qs.  (4.  ZZ)  and  (4.  23),  we  obtain 

AT)  = s.s  SL. 


Comparing  these  results  with  those  of  a plane  normal  shock 
(e.  g.  Ref.  9)t  we  note  first  that  the  plane  shock  strength  ( '\jAYr') 
is  quite  arbitrary  while  fora  cylindrical  shock,  AW 
The  expression  for  shock  thickness  (£q.  4.22)  shows  that 
A.y  ~ 0 (of  ) , although,  combining  , the  expression  (4.  2a) 

agrees  with  that  of  a plane  shock  (cf.  Ref.  9)  within  the  order  of 
magnitude.  The  result  (4.  25)  differs,  however,  from  Lievey's 
result  for  the  diffuse  shock  in  a source  flow  (Red.  4 £q.  4.  9),  in 
which  he  explains  the  discrepancy  as  due  to  some  degree  of  choice 
of  the  definition  of  the  shock  thickness.  Our  result  also  indicates 
that  the  maximum  velocity  gradient  inside  a cylindrical  shock  is 
of  order  <x  , (in  contrast  to  Lievey’s  result:  0(o(  ))  , while 
for  a plane  shock,  the  maximum  gradient  is  (cf.  Ref.  9)  of 
order  which  reduces  to  0(<^  if  AW  0 (cx  '^)  . 

The  expression  (4.  24)  differs  from  the  Prandtl  relation  of  a 

■/j 

plane  shock  by  a term  of  0(o(  ) which  here  agrees  with  Levey's 
resiilt  (Ref.  4,  Eq.  4.  13).  The  differences  between  the  present 
results  of  cylindrical  fl«aw  and  those  of  one -dimensional  plane 
shock  can  perhaps  be  realized  by  visualizing  that  the  viscous 
forces  exerting  on  the  surfaces  r<^  and  of  a fluid  element 
is  indeed  of  quite  different  nature  from  those  exerting  in  plane 
shock  flow,  since  in  the  former  case,  the  normi.!  atress  acting 


I 


(iii) 


(iv) 


On  uuc-fcice  Jr  ivill  have  a component  in  the  radial  direction. 

On  the  subsonic  branch,  w is  a monotonic  decreasing  function 

>/3  2/3 

ot  j . W = / at  point  £ (cf.  Figs.  6 and  4)  where  a 

and  the  velocity  gradient  (0^)  ~~'^£  " ^ 

also  Eq.  2.21a),  which  shows  that  the  solution  curve  of  the 
subsonic  branch  passes  w =/'  (cf.  Fig.  4;  slightly  above  the  curve 
The  thermodynamic  variables  in  this  flow  region  can  be  deduced 
from  Eq.  (4.  5),  (1. 3)  and  (1. 4).  That  is,  in  the  expansion  (4.  3) 
and 


>/3  M ^ ■2/6  M ^ 

Of  />  ($)  r Of  P (^)t  • 

- '/?  o>  2/3  a) 

(?(£)  ? (k)i--- 


(4.  26) 


we  have 


and 


e'h  - -(t-'Kmf }»> 


(4.  27) 


(4.  28) 


where  is  given  in  Eqo.  (4.  14)  and  (4.  16).  The  values  of 

and  8 vs.  ^ is  plotted  in  Fig.  7.  The  supersonic  branch 

starts  with  compression  and  is  then  followed  by  an  expansion 

wave,  while  the  subsonic  branch  expands  continuously,  Eq.  (4.  28) 

states  simply  that  ^ ^ and  9 satisfy  the  isentropic  relation 

up  to  0(<x''^).  This  implies  that  the  entropy  variation,  if  any, 

ty 

across  this  region  must  be  cf  order  at  least  of 
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i 


c 


5.  The  Solution  in  the  Inner  Supeysonic  Region 

The  above  transonic  solution  shows  that  the  supersonic  branch 
flow  approaches  asymptotically  to  the  subsonic  branch  and  for  '’J  <% 
both  branches  have  supersonic  local  speed.  We  shall  proceed  to  find 
the  solution  for  . L.et  us  consider  first  the  continuation  of  the 

subsonic  branch.  For  the  sake  of  convenience,  we  shall  take  the  point 
in  Fig.  6 as  the  boundary  condition  such  that 


w = / 

B 

> % 

and 

(dV 

(5-  1) 


To  obtain  the  solution  for  y^Jg  • <^tie  can  try  the  following  expansion 
vr  = ^ 

y(z)  =.  a y (z,)  -f-  G(  y (X.)  -f  cx  y (z,)  i-  • ‘ ^ (5.2) 

Q{z)  = $ (z)  -f-  <X  $ (z)  f-  ‘ ‘ ‘ 


which  is  to  be  substituted  into  Eqs.  (3.  1)  and  (3.  2).  However,  by  noting 
the  boundary  condition  (5.  1),  a more  convenient  method  to  approximate 
the  solution  can  be  carried  cut  by  letting 
</6 


V = - ^ n a + a a ^ t -h  (^x/z)  fx  i 

df)  « ' \ ^ 


/ ^ 


for 


i = w-/  >o 


(5.3) 


where 


so 


that  the  conditions  (5.  1)  are  satisfied  and,  in  addition 


(5.  3a) 


wiO 


at  'fT  ^ -S  (cf=  ^ 2.2a).  inteeratins  (5. 3)  and 

or  f ' ~ ' 


using  Eq.  (5.  1),  we  obtain 


c 
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, for  i >o 


<5.  4) 


r- 


■ ( 


; = - 


Z a 


»/3 


\jza,a^  -<i< 


^-'hSEb\ 

‘ ^ V.  i£i  ^ 

/ 


This  solution  is  in  good  agreement  with  the  Eq.  (4.  16)  for  ^ • 

It  can  be  seen  at  once  that  as  W ^ ^ (the  maximvim  velocity  at  which 
G^O  1,  or  t y approaches  its  smaller.!  value  , 


tC  C \ 

\^o  -r; 


in  which  the  assigned  value«j  of  ^ ^ ^ have  been  usedo  Thus  we  see 

that  the  flow  supposedly  terminates  itself  at  a distance  of  0(o(  ' ) to 
the  inner  side  of  y-o  ^ beyond  which  there  is  no  solution  to  our  present 
system  of  equations.  To  search  further  for  the  possibilities  whether 
one  still  could  obtain  the  solution  of  physical  reality  for 
would  face  some  rather  dubious  situations.  For  instance,  near  , 

the  density,  temperature  and  pressure  all  become  so  low  that  the  validity 
of  the  equation  of  state  for  perfect  gas  (1.4)  is  questionable.  Besides, 
the  fact  that  the  viscous  stresses  reach  the  magnitude  of  the  fluid  pressuri^ 
near  ^ ^ sets  a likely  limit  as  to  the  applicability  of  the  Navier  - 

Stokes  equation  (1.  1)  and  also  raises  a question  as  to  whether  Burnett's 
higher  viscous  terms  (Ref.  10,  p.  271)  should  be  employed  to  overcome 
the  present  difficulty.  Of  course,  it  would  seem  plausible  to  continue 
our  solution  further  inward  by  assigning  appropriate  values  to  the 
arbitrary  constant  C in  (2.  16).  Nevertheless,  it  is  still  impossible  to 
bring  the  flow  to  ^ = -oo  (r=o)  on  account  of  the  singularity  that  <^u  r ' 
near  r=o  (cf.  1.  3).  To  clarify  these  rather  vague  points  is  beyond  the 
scope  of  this  paper,  although  such  clarification  is  certainly  desirable. 


C 
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6.  Thc.‘  entropy  Variation 


3? 


r 


c 


( 


We  define  S to  be  the  specific  entropy, 

T <Js  Jr  - j-  (6, 1) 

then  the  energy  equation  (1-2)  can  be  written  as 

^7  ti  ^ — J/</  g^rjo  -h  ^ (6.  2) 

where  ^ is  the  viscous  dissipation  function,  which  in  this  case  is, 

p-  = I d/]  . I 3, 

The  above  definition  of  JS  of  a fluid  element  is  clearly  for  an  open  system 
since  the  heat  exchange  by  conduction,  and  hence  a net  flow  of  entropy, 
occurs  with  the  neighboring  elements.  Thus  Eq.  (6.  2)  merely  expresses 
the  energy  balance,  in  terms  of  ^ , of  the  iluid  element  — a system 
not  isolated,  in  the  thermodynamic  sense,  from  its  surroundings.  The 
analysis  of  formulating  the  second  law  of  thermodynamics  for  the  fluid 
flow  case  by  making  the  system  closed  has  been  investigated  in  some 
detail  by  Tolman  and  Fine  (Ref.  11)  and  discussed  later  by  Curtiss  and 
Hirschfelder  (Ref.  12)  from  the  point  of  view  of  statistical  mechanics. 
Their  idea  is  in  essence  to  state  that  the  change  in  the  entropy  of  a 

system  should  be  constituted  of  not  only  the  entropy  carried  into  the 
system  due  to  conduction  of  heat  energy,  equal  to  ^radj)  per 

unit  volume,  but  also  the  net  increase  in  entropy  produced  by  irreversible 
p«rocesses  taking  place  inside  the  system.  Following  Tolman's  notation, 
we  may  thus  write 

^ ^ . ,6.4) 

In  the  present  case,  Eq.  (h.  4),  after  combining  with  (6.  2),  becomes 
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-r^  / ' 


= Y ^ r 


As  we  are  only  interested  in  the  qualitative  features  of  our  later  results, 
we  simplify  these  equations  by  using  the  assumptions 


Y - A^/  ~ ^ ' ^/>  • - constant  • ^ ~Y 


(6,  6) 


Then  the  nondimenaional  form  of  Bqs.  (6.  Z)  anfl  (6.  5)  are  respectively 


4 r (^t-  yf 


dvf 


Of-0  ' \di^J  dTj 


+ 


where 


= % 


and 


'jf-iKd^Jirr  *3  \dyy  dy  J 


(6,  7' 


(6.  7a) 


(6.8) 


Though  the  sign  of  the  terms  ou  the  right  hand  side  of  (6.  7)  is  in  general 
indefinite,  the  value  of  the  l ight  hand  side  terms  of  (6.  8)  is,  however, 
positive  definite.  Therefore  increases  monotonically  along 

the  fluid  flow,  as  predicted  by  the  second  law  for  a closed  system. 
Subtracting  Eq.  (6.8)  from  (6.7),  we  obtain 


a4 

\dyJirr  dy 

which  can  be  integrated  to  yield 


A = d -h  il.  a 

/ 

where  the  constant  of  integration  is  so  chosen  that  both  J and  J.  . 


(6.9) 


tend  to  d as 


, dg  being  arbitrary. 


In  order  to  see  that  d of  the  shock  type  flow  reaches  a maximum 


-19- 


c 


near  w = / , we  substitute  £q.  (2.  3)  into  (6.  7)  and  obtain 


^ ^ 3?  vv  ^ 


y-/  j V dj 


(6.  10) 


This  equation  shows  that  for  y outside  the  transonic  region,  the 


variation  in  ^ is  at  most  of  0(<=<)  . Within  the  transonic  region,  £if 


being  of  0(oT  ) , overwhelms  the  rest  of  the  terms  in  the  bracket  and 
hence  (6. 10)  reduces  to 

y~/  t/»  3 ^ ' 


(6.  11) 


It  then  follows  that 


assumes  its  maximum  value  at  the  point  where 


the  curvature  of  vr  =W(y)  curve  vanishes  ( = o at  point  d in 


z 

dJ 


Fig,  6 and  at  this  point  £f  is  lass  than  zero).  However,  from  (6.9), 

dy^ 

^ d does  not  have  an  extremum  in  the 


the  quantity  J -f-  iL  d(  ^ /tpi 

3 dy  o 

entire  flow  region.  The  above  result  is  very  much  the  same  as  that  of  a 
plane  shock  (e.g.  Ref.  9),  the  solution  of  which  shows  that  the  velocity 
has  a point  of  inflection  at  W:/  where  the  entropy  is  also  maximum. 

Integrating  £q.  (6.  iU)  with  the  aid  of  £qs.  (3.  1)  and  (3.  2)  under 
condition  (6.  6),  we  obtain 


where 


= £ey\^9  i^F)  J =1  Jl^e  -t 

-4  = h (^•/?'> 


(6. 12) 


so  that  J-*-  -Jg  as  F ^ . This  equation  is  actually  the  definition 
of  J usually  given  for  a perfect  gas.  .Substitution  of  the  solution  (3.  12) 
into  (6.  12)  shows  that 

AS  ~ 0(‘'^)  for  y '>  O (/)  (6.13) 

V/ithin  and  around  the  trsneonic  region,  we  substitute  the  solution  (4.  1)- 


-40- 


'sT 


(4.  3)  into  (6.  IZ)  and  simplify  the  expansion,  then  we  obtain 


(6.  14) 


wnere 


( 


(6.  14a) 


and  the  value  ^■=^(x)  is  given  by  Eqs.  (4.  14),  (4.  16)  and  also  plotted  in 

Fig.  6.  Equation  (6.  i4)  is  consistent  with  the  fact  that  ^'  = ^ = const. 

2 

along  the  inviocid  solution  ^ =x  . The  variation  in  4 along  supersonic 
and  subsonic  branches  of  our  solution  follows  directly  from  the  data  shov/n 
in  Fig.  6.  The  remit  is  plotted  in  Fig.  8.  As  y decreases  along  the 
supersonic  branch,  the  entropy  4 first  increases  till  it  reaches  the 

,s/3 

maximum  4^  + 1.  21  or  (^~')  at  point  6 , then  decreases  and 

later  assumes  once  again  the  value  4^  (the  value  of  J at  = « ) at 
point  C where  is  minimum  there.  After  that  4 decreases 

rapidly  with  further  decrease  in  y and  eventually  tend  to  - c>o  as 
the  flow  solution  terminates.  On  the  subsonic  branch,  .4  decreases 
monotically  with  decreasing  f . However,  by  substituting  Eqs.  (6.  14) 
and  (4.  3)  into  (6.  9)«  it  can  easily  be  shown  that  ( 4^)  increases 
monotically  with  decreasing  y and  the  variation  in  (-<(^  ) is  of  order 
. Consequently,  the  result  that  as  y can  be 

explained  by  visualizing  from  £q.  (6.  9)  that  ^ (^9)  decreases 
beyond  all  bounds  as  y-^-y^^  . Physically,  this  probably  implies  that 
the  flow  is  rather  far  from  its  equilibrium  condition  due  to  the  large 
velocity  gradient,  inducing  a rapid  decrease  in  temperature  which  even 
the  important  heat  conduction  in  this  region  can  not  compensate. 


(. 
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7 . The  Case  when  A Proportional  to  Temperature 

In  our  previous  investigation^  ^ and  % were  assumed  to 

be  constant.  In  this  section  the  effects  on  the  solution  due  to  the  variation 
of  y<t  , jji  and  \ , all  assumed  to  be  proportional  to  T , will  be 

calculated  and  compared  with  the  previous  results.  That  is. 


/L  =.  consf.  T , then  Ji  = (3  . {*'.1) 

The  dependence  of  ^ on  T other  than  the  above  relation  can  be  worked 
out  in  a similar  manner.  In  the  present  case  o<  depends  on  TJ  only, 
but  not  on  the  local  temperature  T • Besides,  jul  , ^ and  a all  have 
the  same  dependence  on  T . Hence,  assuming  and  Y still 

to  be  constant,  we  may  take 

i , d . 4.  = , < 


all  constant. 


(7.2) 


Then  the  fundamental  system,  Eqs.  (1.15),  (1.16),  becomes 


+ L \ d (Q  \ _ -_7W  fg  / dif  ^ 


•f- 


^jL*r 


\d6 


(7.  3) 


and 


*v  , 9 

T *j:rr*w 


'Oi-i)  dv  j 


Yh 

zCY-0 


(7.4) 


L dy  (r-Od'O^l)  <^y 

In  the  outer  region,  we  substitute  the  same  expansion  Eq.  (3.  4) 
into  Eqs.  (7.  3)  and  (7.4),  then  carry  out  the  calculation  in  a way  similar 
to  that  described  in  % 3.  The  result  is  found  to  be 

«(.i)  = i + I’l] 


(7.5) 


where 


A - izE  - 


-yo.uJ,  . „.sa, 


2 2 r ' ' V 

and  y(^)  has  the  same  expression  as  Eq.  (3.  12b)  up  to  terms  of 
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significant  order.  Comparing  this  solution  with  £q.  (3.  IZ),  one  can  see 


that  for  either  H-  ^ » or  0<y/r  < /-0(o<^^)  the  correction 

to  the  value  of  w at  ^ due  to  the  variation  in  ^ is  at  most  of  order 
. A correction  of  the  same  order  also  applies  to  6 . Thus 
w^e  may  conclude  that  the  effect  cf  varying  ^ is  rather  unimportant  in 
this  range. 

In  the  transonic  region,  we  can  again  try  to  obtain  the  solution 
of  Eqs.  (7.  3),  (7.  4)  in  the  form  of  Eqs.  (4.  1)-  (4.  3).  This  substitution 
leads  to  the  following  equations: 

f/ - ^ ^ ^ ^ (7.6; 

y/ 


'/j. 


0- 


- i-6Q-3/f)  g 

/ -« 


,(f'7 


0)  .a 


1H  / 


1 


and 


9(S)  = /-  -2  'I  + 


rOfi)  J , 


0(c() 


(1.1) 


(7o8) 


Putting  first  in  these  equations  and  then  comparring  them  with 

the  corresponding  equations  (4.  5)-  (4.  7)  for  the  case  of  constant  jl  , we 
see  that  the  effect  of  varying  /t  is  the  introduction  of  an  extra  term 
to  the  second  order  equation,  namely,  the  first  term  on  the  right  hand 
side  of  Eq.  (7.  7).  This  again  implies  that  to  account  for  the  effect  of 
variation  in  yH  , the  values  of  tv  and  0 should  be  i;ubjected  to  a 
correction  term  of  order  0(<*  ) . It  can  also  be  shouna  that  in  the 

inner  supersonic  region,  its  effect  is  adso  of  0 ) . For  the  case 
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where  the  dependence  of  ^ on  T is  different  from  Eq.  {7.  1)  (for 
instance  JU.  - C t'^  , ^ > o ),  the  effect  should  qualitatively  remain 

the  same. 


( 


c 
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